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ABSTRACT 

We present results on the consistency of standard cold dark matter (CDM) models 
with both the COBE normalization and the data from the 4th channel of the 9-point South 
Pole scan. We find that CDM models are consistent with both experiments, a conclusion 
which is at odds with some other analyses. This is partly due to our careful treatment 
of the temperature autocorrelation function, but also derives from a realization that the 
statistical conclusions depend strongly on assumptions about the prior distribution for the 
normalization Q. 

Subject headings: cosmic microwave background — cosmology: theory 
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1. Introduction 



There has been some confusion in the literature concerning the implications for cold 
dark matter (CDM) models of the results of the COBE (Smoot et al. 1992) and SP91 (Gaier 
et al. 1992) experiments. In particular, there have been contradictory claims (e.g. Dodelson 
& Jubas 1993, Gorski et al. 1993) about the fundamental question: how incompatible are 
standard CDM models with these anisotropy experiments? Our intention is to provide a 
definitive answer to this question, giving the details of our procedure. 

We analyze in detail the data from the 4th (highest frequency) channel of the SP91 9- 
point scan (Gaier et al. 1992). This data set seems most likely to represent actual cosmic 
microwave background (CMB) fluctuations, and is the data set previously analyzed by 
Gorski (1992), Dodelson & Jubas (1993) and Gorski et al. (1993). We wish to stress three 
points about our analysis: 

(1) We use the theoretical two-point autocorrelation function predicted for the SP91 
measurement strategy by standard CDM models (fi = 1, h = 0.5, n = 1, 0.01 < fi B < 
0.10, A = 0, T/S = 0). We notice that points whose angular separation is near the peak- 
to-peak chopping angle have strong negative correlations. This result, explained below, is 
entirely due to chopping, and is completely independent of the underlying cosmological 
model. We compare the CDM autocorrelation function with the autocorrelation function 
assumed by Gaier et al. (1992) in their analysis. 

(2) We quote results for an experiment-independent parameter, which we take to be 
Q = Qrms-ps = <<5|ms>°' 5 (Smoot et al. 1992, Wright et al. 1993), the r.m.s. value of 
the quadrupole moment averaged over an ensemble of universes. In contrast, the usually 
quoted values of Co as a function of C (see below) must be disentangled from the observing 
strategy to find information on the underlying spectrum of fluctuations. Deconvolving 
theory from experiment is essential in obtaining a meaningful comparison of experiments 
with different strategies measuring different parts of the sky. We remind the reader that 
the value of Q measured by COBE (Smoot et al. 1992, Wright et al. 1993) from the overall 
normalization of the experimental two-point correlation function (assuming CDM with 
n = 1 and no gravity wave contribution) is Q = 17 ± 5 /iK. The quoted error includes the 
effects of cosmic and sample variance (e.g. White et al. 1993, Scott et al. 1993). 

(3) We perform a maximum likelihood analysis, but investigate sensitivity to the "prior 
distribution"; that is, should we assume that, a priori, equal intervals in Q are equally 
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likely (the usual assumption), or equal intervals in Q 2 (which is simply proportional to 
the amplitude of the power spectrum) , or some other choice? Different choices turn out to 
change the upper limit on Q from SP91 (at various confidence levels) by tens of per cent. 
Note that the choice of prior becomes less important in the limit that the data have great 
"discriminating power". At present, sensitivity to the prior can be used as a measure of 
how constraining the 9-point scan is. 

2. The Temperature Autocorrelation Function 

We begin with the theoretical two-point autocorrelation function for CDM models. Let 
T(6>, 4>) denote the temperature difference that the experiment assigns to a point on the 
sky: 

T(9, <P)=QJ2 a ^ Wtm Y tm (0, (j>) , (1) 
£m 

where the W^ m 's represent the window function of the experiment, and the a^ m 's are 
random variables whose distribution must be specified by a specific cosmological model. 
In general, rotational invariance implies that 

(at m a\, m ,) = Ce8ee> b~ mm ' , (2) 

where the angular brackets denote an ensemble average over the probability distribution 
for the a^ m 's, and Cg is normalized so that = 4-7r/5. For a pure Sachs- Wolfe, n = 1 
spectrum, C~[ l oc £(£ + 1). We compute the C/s for CDM models using power spectra 
provided by Sugiyama (e.g. Sugiyama & Gouda 1992), which are essentially identical with 
those computed by Bond & Efstathiou (e.g. 1987). 

The theoretical two-point autocorrelation function is given by 
C t h,y = Q~ 2 (TiOi^dTiOj,^)) 

= J2°e\ W ^\ 2 Y em (0i, <t>i) YLiOj, <t>j) , (3) 

where we have divided by Q 2 to make C t h dimensionless; there is then no Q dependence 
at all in C t h- Note that C t h is calculated using an underlying theory (e.g. CDM) together 
with the experimental window function, and does not depend on any experimental data. 
Thus C t h is not the same thing as the experimental correlation function measured by the 
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COBE team, which is computed by taking a sky average of data without any theoretical 
input. 

All the points in the 9-point scan of SP91 have the same zenith angle but different 
azimuth. In (3) we can write 9i = 9 = 27°75 and (pi = Acpxi with A<psm9 = 2°1 being the 
angular spacing of points on the sky. Furthermore the chopping is done in the direction, 
which implies that 

W tm = 2H (ma) exp[-(£ + y 2 ) V/2] (4) 

(Bond et al. 1991, Dodelson & Jubas 1993, White et al. 1993), where H is the Struve 
function (Gradshteyn & Ryzhik 1980), asin6> = 1°5 is half the peak-to-peak chopping 
angle on the sky, and a = 0.425 x 1°4 is the Gaussian beam width of the antenna. Thus, 
equation (3) becomes 

oo I 

Cth.ij = J2 Ce ex Ph(^+ V2)V 2 ] J2 4 #oV«) IW#,0)| 2 cos(m^) , (5) 

£=2 m=-£ 

where 0^ = Ac/) x \i — j\. We find that C t h is sensitive to the value of a, and expect that 
departures from a Gaussian beam profile may also significantly affect the results. 

In Figure 1, we show C t h as a function of 4>sm9, the separation angle on the sky, for 
fi B = 0.03. We see that it becomes strongly negative for <psin9 near 3°, corresponding 
to the peak-to-peak chopping angle on the sky, 2a sin 9. This is easy to understand: the 
region halfway between the two negatively correlated points contributes positively to one 
point and negatively to the other, which will clearly produce an anticorrelation. 

Results of small-scale CMB experiments are often quoted assuming an underlying 
Gaussian autocorrelation function (GACF); this corresponds to taking 

Q 2 C i = 2ir9 2 c C exp[-(£+ y 2 ) 2 9 2 c /2], (6) 

where 9 C and Cq are parameters. In Figure 1, we show the corresponding Cth as a function 
of <£sin0 for 9 C = 1°2 and C /Q 2 = 7.9. The choice of 9 C = 1°2 is the "best-fit" value of 
Gaier et al. (1992), which matches the peak of equation (6) to the peak of the SP91 window 
function in multipole space. The choice of Cq/Q 2 = 7.9 results in C t h(0) = 4.9, the same 
value as in the CDM model with fi B = 0.03. Notice that for these particular choices of 9 C 
and C , the anti-correlation at the nearest-neighbor separation of A</>sin6> = 2°1 is well 
reproduced. This agreement would not hold for general values of 9 C and Cq. 
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The signal-to-noise of the experiment is currently such that experimental errors con- 
tribute large diagonal entries to the correlation matrix [see equation (8) below], reducing 
the effect of off-diagonal entries on the fit. If the signal-to-noise were to improve signif- 
icantly, then the off-diagonal entries would play a larger role. As an example, with the 
current data we find that a diagonal Cth,ij produces a very similar limit to that obtained 
in Gaier et al. (1992), showing that the assumed off-diagonal terms are playing little role 
in the analysis. This can also be seen in Table 1, where the upper limits on Q for a given 
choice of prior distribution (see below) simply scale like C th ' (0) as fi B is varied. 



3. Maximum Likelihood Analysis 

We now turn to the limits which can be placed on Q. In accord with standard CDM 
models we consider underlying cosmological fluctuations with a Gaussian probability dis- 
tribution (not to be confused with fluctuations having a Gaussian power spectrum, or 
GACF). We assume that the experimental errors cr; for each Tj are uncorrelated and 
Gaussian distributed. Then the unnormalized likelihood function for Q is given by 

£(Q) « -/^expf-VsT^- 1 )^] , (7) 

where the matrix K is 

Kij = Q 2 C thjij + of Sij . (8) 

Equation (7) assumes that the temperatures have no systematic errors. However, SP91 
has a possible, unknown, systematic offset and linear gradient, and so their best-fit values 
have been removed from the measured temperatures. [In fact, as a consequence of the 
data binning procedure, the published values of Tj (Gaier et al. 1992) have a small residual 
gradient which should be removed (Gaier 1993).] To project out a weighted offset and 
linear gradient from the Tj, we write T a = R a iTi where R is a 7 x 9 matrix (a = 1, . . . , 7 
and i = 1, . . . , 9). The seven rows of R must be linearly independent, and each row must 
be orthogonal to each of two nine-component vectors: one whose elements are 1/of , and 
one whose elements are (i — 5)/of . The choice of R is otherwise arbitrary and does not 
affect the results. Then, instead of equation (7), the unnormalized likelihood function is 
given by 

C(Q) oc ^ = ^exp[-V2T a (M- 1 ) a6 f & ] , (9) 
V det M 
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where M ab = R ai KijRj b . Note that the procedure used by Bond et al. (1991) is equivalent 
to removing an unweighted offset and linear gradient from the data. 

We can now use equation (9) to set confidence levels on Q. We must first choose a 
prior distribution for Q; that is, we must decide whether equal intervals of Q, Q 2 , or some 
other monotonic function f(Q), are equally likely a priori. Given a choice of f(Q), the 
upper limit on Q at a confidence level of c is the solution Q max to the equation 

_ J Q — C(Q)df(Q) 
C J-C(Q)df(Q) ■ UUJ 
Due to the large measurement errors in the data set, assuming f(Q) = logQ, as recom- 
mended for an overall scale factor, can be problematic (see e.g. Readhead et al. 1989). One 
is thus led to consider prior distributions which are uniform in the scaling variable. The 
usual choice (Bond et al. 1991, Dodelson & Jubas 1993, Gorski et al. 1993) is f(Q) = Q 
(the "bias" parameter b p is proportional to so that dQ = db p /b 2 p ). However, there is 
no compelling reason to assume a prior distribution that is uniform in Q. It is just as nat- 
ural, for example, to assume that the prior distribution is uniform in the power spectrum 
normalization Q 2 , or to make use of the COBE measurement of Q. The choice of f(Q) is 
important because it strongly affects the resulting values of <5 ma x- Note that for "good" 
data, the choice of the prior distribution should make very little difference, so "prior de- 
pendence" gives us a handle on the constraining power of the data. This is illustrated in 
Figure 2, where we show the cumulative likelihood (c vs. Qmax) for four different choices 
of fi B and for f(Q) = Q and f(Q) = Q 2 ■ Clearly the effect of the prior distribution is 
significant. This sensitivity may be understood as being due to limited sky coverage (see 
Scott et al. 1993). 

In Figure 2 we also show a band corresponding to the COBE measurement of Q = 
17 ± 5/uK (la errors). We feel that this plot makes it clear that CDM models cannot 
be ruled out at high levels of confidence from a combination of the COBE and SP91 9- 
point data. This is also shown in Table 1, where we list the upper limits on Q at the 95% 
confidence level for the eight combinations of Q B and choice of prior distribution f(Q) = Q 
and f(Q) = Q 2 . 

4. Conclusions 

A reanalysis of the SP91 data from the 9-point scan shows no conflict with standard 
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CDM models for any value of fl B and for Q within the one-sigma range specified by 
the COBE results, Q = 17 ± 5 juK, in contradiction with the conclusions of some earlier 
analyses. We find that the window function for SP91 is sensitive to the beam size assumed 
and suspect that deviations from a Gaussian beam pattern may significantly affect the 
results. 

We are extremely grateful to Naoshi Sugiyama for providing us with CDM radiation 
power spectra. We also thank Josh Gundersen, Phil Lubin, Peter Meinhold, and especially 
Todd Gaier for helpful discussions. E.B. acknowledges the support of an NSF graduate 
fellowship. This work was supported in part by NSF Grant Nos. PHY-91-16964 and 
AST-91-20005. 
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16 


24 


0.06 


5.4 


15 


22 


0.10 


6.2 


14 


21 



Table 1: The 95% confidence level upper limits, Q max (a*K), from the 9-point scan for 
CDM with a range of Vt B and for prior distributions dQ and d(Q 2 ). Also listed is the value 
of Cth(0) for each Q B . 
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Figure Captions 



Fig. 1. The theoretical temperature autocorrelation function, Cth, as a function of scan 
angle on the sky, <psm9 (<p is the azimuthal angle and 9 the zenith angle). The solid line 
is the autocorrelation function assuming the CDM power spectrum for Q B = 0.03 and 
h = 0.5. The dashed line is the traditional Gaussian autocorrelation function (GACF) 
approximation, with a correlation angle 9 C = 1°2, processed through the experimental 
observing strategy and matched to the CDM prediction at </> = 0. For comparison, the 
corresponding power spectra, l[l + 1)CV, for both theories are shown in the inset. 

Fig. 2. The cumulative likelihood as a function of power spectrum normalization, Q (/LtK), 
for CDM models with Q B = 0.01,0.03,0.06,0.10 (fi B increases to the left). Curves are 
labelled by the assumed prior distribution. The hatched band is the ±lcr allowed range of 
Q from COBE. 
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